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Abstract 

Wc give new sufficient conditions for a sequence of polynomials to have only real zeros 
based on the method of interlacing zeros. As applications we derive several well-known facts, 
including the reality of zeros of orthogonal polynomials, matching polynomials, Narayana 
polynomials and Eulcrian polynomials. We also settle certain conjectures of Stahl on genus 
polynomials by proving them for certain classes of graphs, while showing that they are false 
in general. 
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1 Introduction 

Polynomials with only real zeros arise often in combinatorics and other branches of math- 
ematics (see [SI [TUl EEJ EOj ESI HH SSI ESI ESI El])- Since our interest is combinatorial, we are 
mainly concerned with polynomials whose coefficients are positive or alternating in sign. The 
real zeros of such polynomials are either all negative or all positive. A polynomial with coeffi- 
cients alternating in sign can be converted into a polynomial with positive coefficients. So we 
may concentrate our attention on polynomials with positive coefficients. We were led to the 
study of such polynomials because of their implications on unimodality and log-concavity. 

Let ao, a±, . . . , a n be a sequence of positive real numbers. The sequence is said to be unimodal 
if there exists an index < m < n, called the mode of the sequence, such that oq < • • • < 
a-m-i < a-m > flm+i > • • • > a n - The sequence is said to be log-concave if ai-\a. l+ i < a\ for 
% = 1, . . . , n— 1. Clearly, log-concavity implies unimodality. Unimodal and log-concave sequences 
occur naturally in combinatorics, analysis, algebra, geometry, probability and statistics. The 
reader is referred to Stanley [44] and Brenti p] for surveys and ^4] [TSl [57l [58l [59l [60l IHT] [62] 
for recent progress on this subject. 

One classical approach to unimodality and log-concavity of a finite sequence is to use New- 
ton's inequality: if the polynomial Y17=o a ^ x% with positive coefficients has only real zeros, then 
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for 1 < i < n — 1, and the sequence ao,a±, . . . ,a n is therefore unimodal and log-concave (see 
Hardy, Littlewood and Polya (3TJ p. 104]). Such a sequence of positive numbers whose generating 
function has only real zeros is called a Polya frequency sequence in the theory of total positivity. 
See Karlin [M] for a standard reference on total positivity and Brenti [10^ [T4"l [T5] for applications 
of total positivity to unimodality and log-concavity problems. It often occurs that unimodality 
of a sequence is known, yet to determine the exact number and location of modes is a much 
more difficult task. The case for Polya frequency sequences is somewhat different. Darroch [19] 
showed that if the polynomial P{x) = Y17=o aiX% w ith positive coefficients has only real zeros, 
then the unimodal sequence oq, a±, . . . , a n has at most two modes and each mode m satisfies 



< m < 
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Our main concern here is sequences of polynomials with only real zeros. Such polynomial 
sequences occurring in combinatorics often satisfy certain recurrence relations. For example, the 
sequence of classical orthogonal polynomials p n (x) satisfies a three-term recurrence relation 

p n (x) = (a n x + b n )p n -i{x) - c n p n - 2 {x) (1.1) 

with p„i(x) = and po(x) = 1, where a n ,c n > and b n G M (see [53j Theorem 3.2.1]). A 
standard result in the theory of orthogonal polynomials is that for each n > 1, the zeros of p n ( x ) 
are real, simple, and separate those of p n+ i(x). In this paper, we develop methods to provide 
a unified approach to the reality of zeros of polynomial sequences satisfying certain recurrence 
relations. 

Following Wagner [56] , a real polynomial is said to be standard if either it is identically zero or 
its leading coefficient is positive. Let RZ denote the set of real polynomials with only real zeros. 
In particular, let PF consist of those polynomials in RZ whose coefficients are nonnegative. In 
other words, PF is the set of polynomials whose coefficients form a Polya frequency sequence. 
Clearly, all zeros of a polynomial in PF are real and nonpositive. For convenience, let G PF. 

Suppose that f,g G RZ. Let {r^} and {sj} be all zeros of / and g in nonincreasing order 
respectively. We say that g alternates left of f (g alternates / for short) if degf = degg = n 
and 

s n <r n < s n _i < • • ■ < s 2 < r 2 < si < r%. (1.2) 
We say that g interlaces f if deg / = deg g + 1 = n and 

r n < Sn-l < ■ ■ ■ < S2 <r 2 < si < n. (1.3) 

Let g ^ / denote "either g alternates / or g interlaces /". If no equality sign occurs in (jl.2p 
(resp. (|1.3j) ). then we say that g strictly alternates / (resp. g strictly interlaces /). Let g -< / 
denote "either g strictly alternates / or g strictly interlaces /". For notational convenience, let 
a < bx + c for any real constants a, b, c and / ^ 0, ^ / for any real polynomial / with only 
real zeros. 

Let {P n (x)} n >o be a sequence of standard polynomials. We say that {P n (x)} is a Sturm 
sequence if degP n = n and P n -i(r)P n+ \{r) < whenever P n (r) = and n > 1. It is well 
known that {P n (x)} is a Sturm sequence if and only if, for n > 1, P n G RZ and P n strictly 
interlaces P n +i- We say that {P n (x)} is a generalized Sturm sequence if P n G RZ and Pq -< P\ ■< 
■ ■ ■ H P n -i ^ P n -< ■ ■ ■ . For example, if P is a standard polynomial with only real zeros and 
deg-P = n, then P", p( n_1 ) 5 . . . } p form a generalized Sturm sequence by Rolle's theorem. 
In particular, if the zeros of P are distinct, then p(«) j p(«-i) ) ...,P',P form a Sturm sequence. 

There are various methods for showing that polynomials have only real zeros. For example, 
Wang and Yeh |61| Theorem 1] established the following result which has been proved to be an 
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extremely useful tool. In fact, it provides a unified approach to unimodality and log-concavity 
of many well-known sequences in combinatorics. See [61] for details. 

Theorem 1.1. Let f and g be real polynomials whose leading coefficients have the same sign. 
Suppose that f,g£ RZ and g ■< f . If ad < be, then (ax + b)f(x) + (cx + d)g(x) G RZ. 

An immediate consequence of Theorem 11.11 is the following corollary ( |61| Corollary 1]). 

Corollary 1.2. Suppose that f,g G PF and g interlaces f. If ad > be, then (ax + b)f(x) + 
x(cx + d)g(x) G RZ. 

Theorem 11.11 and Corollary 11.21 provide the inductive basis for showing the reality of zeros of 
polynomial sequences {P n (x)} satisfying certain recurrence relations P n (x) = a n (x)P n _i(x) + 
b n (x)P! a _ l (x). For example, let S(n, k) be the Stirling number of the second kind and let S n (x) = 
Ylk=o S( n , k)x k denote the associated generating function. It is well known that S(n, k) = 
S(n — 1, k — 1) + kS(n — l,k), S(0, 0) = 1 and S(n, 0) = for n > 1, which is equivalent to 

S n (x) = xS n -i(x) + xS'n^x) (1.4) 

with Sq(x) = 1. From Corollary 11.21 (or Theorem I l.ip and by induction, it follows immediately 
that S n (x) has only real zeros for each n > 1, a result originally due to Harper [32]. Similarly, 
the classical Eulerian polynomial A n (x), which satisfies the recurrence relation 

A n (x) = nxA n _i(x) + x(l — x)A' n _ 1 (x) (1-5) 

with Aq(x) = 1 ( see Comtet \17\ Exercise VII. 3] for instance), has only real zeros for n > 1. 

However, Theorem 11.11 and Corollary 11.21 cannot be used to prove the fact that {S n (x)} and 
{A n (x)} form generalized Sturm sequences respectively. Neither are applicable to polynomial se- 
quences {P n (x)} satisfying recurrence relations similar to P n (x) = a n (x)P n -\(x) + b n (x)P n -2(x) , 
nor even to the orthogonal polynomials. So a natural problem arises: Given F(x) = a(x)f(x) + 
b(x)g(x), under what conditions g ^ / necessarily implies / ^ Fl This leads to the following 
two more general problems. 

Problem 1.3. Let F(x) = a(x)f(x) + Y^j=i bj( x )dj( x )- Suppose that f,gj G RZ and gj H / for 
all j. Under what conditions F G RZ and f ^ F? 

Problem 1.4. Let F(x) = a(x)f(x) + b(x)g(x) and G(x) = c(x)f(x) + d(x)g(x). Suppose that 
f,g£ RZ and g ^ / . Under what conditions F,G G RZ and G ^ F? 

The organization of this paper is as follows. In Section 2, we present various results concern- 
ing Problem II. 31 and II .41 Among other things, we give a simple proof of Theorem ll.il In Section 
3, we apply these results to derive several well-known facts and solve certain open problems and 
conjectures in a unified manner. 



2 Main results 



Let sgn denote the sign function defined on M by 

'+1 if cc > 0, 
sgn (x) = < if x = 0, 
-1 ifrc<0. 

Let f(x) be a real function. Denote sgn/(+oo) = +1 (resp. 
for sufficiently large x. The meaning of sgn/(— oo) is similar. 



■1) if sgn/(x) = +1 (resp. — 1) 
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Theorem 2.1. Let F,f,g be three real polynomials satisfying the following conditions. 

(a) F(x) = a(x)f(x) + b(x)g(x), where a(x), b(x) are two real polynomials, such that degF = 
deg / or deg / + 1. 

(b) /,<?GRZ andg<f. 

(c) F and g have leading coefficients of the same sign. 

Suppose that b(r) < whenever f(r) = 0. Then F G RZ and f ^ F. In particular, if g < f and 
b(r) < whenever f{r) = 0, then f ~< F. 

Proof. Without loss of generality, we may assume that / and g have no common zeros, which 
implies that / has only simple zeros. In other words, it suffices to consider the case g -< / '. We 
may also assume that F, g are standard and that dega(x) < l,deg6(x) < 2. If b(x) = 0, then 
the result is trivial. So let b(x) ^= 0. 

First consider the case b(r) < whenever f(r) = 0. Let deg/ = n and r n <•••<. r\ 
be all zeros of /. Then sgng(rk) = (— l) k ~ 1 since g is standard and g -< /. It follows that 
sgnF(rfc) = (— l) k . Since F is standard, we have sgnF(+oo) = +1 and sgnF(— oo) = (— l) n+1 
provided degF = n + 1. By Weierstrass Intermediate Value Theorem, F(x) has one zero in 
each of n intervals (r n , r n -±), . . . , (7*2, r%), (r\, +00) and has an additional zero in the interval 
(— 00, r n ) provided degF = n + 1. Hence F has degF real zeros and f < F. 

For the general case, define bj(x) = b(x) — 1/j and Fj(x) = a(x)f(x) + bj(x)g(x). Let j 
be sufficiently large. Then bj{rk) < for each zero of / since the number of zeros of b(x) 
is finite, and so Fj G RZ and / -< Fj by the above discussion. It is clear that degFj = degF 

since Fj = F — g/j. Let r\ > r2 > r$ > • • • and t± > 4 > ^3 > ' ' ' be all zeros of / 

(i) (i) (i) 

and Fj respectively. Then t± > r\ > t% > r% > > r% > ■ ■ ■ . Note that the zeros of a 

polynomial are continuous functions of the coefficients of the polynomial (see |18j for instance). 

In particular the limit of a sequence of RZ polynomials is still a RZ polynomial. Hence F G RZ. 

Moreover, let t\ > ti > £3 > • • • be all zeros of F. Then t\ > r\ > t<i > r<i > t$ > r% > ■ ■ ■ by 

continuity. Thus / ^ F and the proof is complete. □ 

Generally speaking, the conditions (a), (b) and (c) in Theorem 12. II can be satisfied naturally. 
It remains to examine the sign of b(x) for the zeros of /. Sometimes this task is trivial, for 
example, when b(x) = —(bo + b\x) 2 . On the other hand, if coefficients of / are nonnegative 
(resp. alternating in sign), then it suffices to consider the sign of b(x) for x < (resp. x > 0). 
As an example, we give a short and simple proof of Haglund [301 Lemma 3.6]. We also refer the 
reader to Theorem 2.4.2-2.4.6 in Brenti [10] for which Theorem 12 . 1 1 can be used to give unified 
proofs. 

Corollary 2.2 ( \30\ Lemma 3.6]). Let f and g be two real polynomials with positive leading 
coefficients a and (3 respectively. Suppose that the following conditions are satisfied. 

(a) f,g G RZ and g interlaces f. 

(b) F(x) = (ax + b)f(x) + x(x + d)g(x) where a,b,d 6i with d > 0, d > b/a, and either a > 
or a < —(3 j a. 

(c) All zeros of f are nonpositive if a > and nonnegative if a < — (3 /a. 

Then F G RZ. In addition, if each zero r of f satisfies —d < r < 0, then f interlaces F. 
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Proof. Suppose that a > 0. Then / G PF, and so g G PF. It follows from Corollary 11.21 that 
F G RZ since a ■ d > b ■ 1 . 

Now suppose that a < —[3 /a. Then / and g have coefficients alternating in sign and F has 
negative leading coefficient. Let deg/ = n. Then degg = n — 1 and deg-F = n + 1. Define 
^(x) = (_i)"/(_ x ) 59l ( x ) = an d = (_i)ni?(_ x ). Then f 1 ,g 1 G PF and 

5(1 interlaces /1. Note that Fi(x) = (—ax + b)fi(x) + x(— x + d)gi(x) and (—a) ■ d > b ■ (—1). 
Hence F x G RZ by Corollary O and so F G RZ. 

Finally, if —d < r < whenever /(r) = 0, then /, g G PF and r(r + d) < 0. Thus / interlaces 
F by Theorem EQ □ 

The following theorem is a generalization of Theorem 12.11 and gives a solution to Problem ll.31 
It can be proved by the same technique used in the proof of Theorem 12.11 So we omit its proof 
for brevity. 

Theorem 2.3. Let F, f,g±, . . . ,gj. be real polynomials satisfying the following conditions. 

(a) F(x) = a(x)f(x) + b\(x)gi(x) + ••• + bk(x)gk{x), where a(x), b\(x), . . . , bk(x) are real 
polynomials, such that degF = deg/ or deg / + 1. 

(b ) f, 9j G RZ and gj < f for each j . 

(c) F and gi,--.,gk have leading coefficients of the same sign. 

Suppose that bj(r) < for each j and each zero r of f . Then F G RZ and f ^ F. In particular, 
if for each zero r of f , there is an index j such that gj -< / and bj(r) < 0, then f -< F. 

Corollary 2.4. Let {P n (x)} be a sequence of standard polynomials and satisfy the recurrence 
relation 

P n (x) = a n (x)P n -i(x) +b n (x)P' n _ 1 {x) + c n {x)P n _ 2 {x), 

where a n (x),b n (x),c n (x) are real polynomials such that degP n = degP n ,_i or deg-P n ,_i + 1. 
Suppose that for each n, coefficients of P n (x) are nonnegative (resp. alternating in sign). If 
b n (x) < and c n (x) < whenever x < (resp. x > 0), then {P n (x)} forms a generalized Sturm 
sequence. In particular, if for each n, degP n = n and either b n (x) < or c n (x) < whenever 
x < (resp. x > 0), then {P n (x)} forms a Sturm sequence. 

Corollary 12.41 provides a unified approach to the reality of zeros of certain well-known poly- 
nomials, including the orthogonal polynomials p n (x) satisfying (II. ip . the Stirling polynomials 
S n (x) of the second kind satisfying (jl.4p and the Eulerian polynomials A n (x) satisfying (jl.5p . 
We will give some more applications of Corollary 12.41 in the next section. 

Lemma 2.5. Let G(x) = c(x)f(x) + d(x)g(x) where G,f,g are standard and c(x),d(x) are real 
polynomials. Suppose that f,g£ RZ and g < f ■ Then the following hold. 

(i) If deg G < degg + 1 and c(s) > whenever g(s) = 0, then G G RZ and g -< G. 

(ii) If deg G < deg / and d(r) > whenever f(r) = 0, then G G RZ and G < f . 

The statements also hold if all instances of -< and > are replaced by ^ and > respectively. 

Proof, (i) Let deg g = m and s m < s m _i < • • • < s 2 < si be all zeros of g. Then sgnG(sfc) = 
(— l) fe since c(sk) > 0. Also, sgnG(+oo) = +1, and sgnG(— 00) = (— l) m+1 provided degG = 
m + 1. Hence G has one zero in each of in intervals ), . . . , (s 2 , -si), (si, +00) and has 

an additional zero in the interval (—00, s m ) provided degG = m + 1. Thus G G RZ and g -< G. 
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(ii) Let deg/ = n and r n < r n _i < • • • < T2 < r\ be all zeros of /. Then sgn G(rfc) = 
(— and sgnG(— oo) = (—1)™ provided degG = n. It follows that G has one zero in each 
of n — 1 intervals (r n ,r n _i), (r^ri) and has an additional zero in the interval (—00, r n ) 
provided degG = n. Thus G G RZ and G -< f. 

The remaining part of the lemma follows by a continuity argument. □ 

The following theorem gives a solution to Problem 11.41 

Theorem 2.6. Let f,g,F,G be four standard real polynomials satisfying the following condi- 
tions. 

(a) F(x) = a(x)f(x) + b(x)g(x) and G(x) = c(x)f(x) + d(x)g(x) where a(x), b(x), c(x), d(x) 
are real polynomials such that deg F = deg G or deg G + 1 . 

(b) /,<?GRZ andg<f. 

(c) A(x) := a(x)d(x) — b(x)c(x) > whenever G(x) = 0. 

Suppose that either c(x) is a positive constant and deg G < degg+l or d(x) is a positive constant 
and deg G < deg /. Then F,G £ RZ and G ^ F. In particular, if g -< / and A(x) > whenever 
G{x) = 0, i/ten G -< F. 

Proof. Suppose that c(x) is a positive constant. Then by Lemma 12.51 (i), g ^ / (resp. g < f) 
implies that G G RZ and g ^ G (resp. g ~< G). By Condition (a), we have oF = aG+ (be — ad)g. 
If A(x) > (resp. A(x) > 0) whenever G(x) = 0, then by Theorem EH F G RZ and G H F 
(resp. G -< F). 

Suppose that d(x) is a positive constant. Then by Lemma [2. 51 (h). g ^ f (resp. g < f) implies 
that G G RZ and G ^ / (resp. G ^ /). By Condition (a), we have dF = (ad - bc)f + bG. If 
A(x) > (resp. A(x) > 0) whenever G(x) = 0, then by Lemma [231 (i), F G RZ and G < F 
(resp. G -< F). □ 

When both c and d are constants, Theorem 12.61 is particularly interesting and useful as we 
shall see. 

Corollary 2.7. Let f,g,F,G be standard real polynomials and satisfy the following conditions. 

(a) F(x) = a(x)f(x)+b(x)g(x) andG(x) = cf(x)+dg(x) where a(x),b(x) G R[x] andc,d G M. 

(b ) deg F = deg G or deg G + 1 . 

(c) f,g£ RZ and^f ^ /. 

Suppose that da(x) > cb(x) whenever G(x) = 0. Then F,G G RZ and G < F. In particular, if 
g < f and da(x) > cb(x) whenever G(x) = 0, then G < F . 

It is well known that if g ■< f, then cf + dg G RZ for arbitrary real numbers c and d, which 
is an immediate consequence of Theorem 11.11 and can also be obtained by setting a(x) = c and 
b(x) = d in Corollary 12.71 Using Corollary 12.71 we can obtain more precise results. The following 
are special cases of Corollary 12 . 71 when both a and b are constants. Some of them have appeared 
in the literature, for example, Wagner |551 156j . 

Corollary 2.8. Let a,b,c,d > 0. Suppose that f,g G RZ are standard and g H /. Then the 
following statements hold. 

(i) If ad> be, then cf + dg ^ af + bg. In particular, g ^ af + bg ^ / . 
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(ii) If af — bg is standard, then cf + dg ■< af — bg. In particular, f,g ^ af — bg. 

(Hi) If —af + bg is standard, then —af + bg^cf + dg, and in particular, —af + bg ■< f,g. 

Similar results hold when g -< /. 

Using Corollary 12.71 we can give a short and simple proof of Theorem II. 11 In fact, we can 
prove the following result more precise than Theorem 11.11 

Corollary 2.9. let f,g £ RZ have leading coefficients of the same sign and g ^ /. Define 
F = (ax + b)f + (cx + d)g and G = af + eg where a, b,c,d£ R. // ad < be, then F, G £ RZ and 
G^F. 

Proof. We have F = xG + H where H = bf + dg. Clearly, G,H £ RZ. If G = 0, then the 
statement is trivial, so let G ^ 0. Without loss of generality, we may assume that both / and g 
are monic. We may also assume that G is standard (otherwise replaced G by —G). Note that 
(ax + b)c — (cx + d)a = be — ad > 0. To prove the result by means of Corollary 12.71 it suffices 
to prove that F is standard. 

Actually, if a = or degG = deg/, then F is obviously standard. Assume now that a ^ 
and degG < deg/. Then deg/ = degg and a = — c. Let / = YYl=i( x ~ r i) an d 9 = ]Xi=i( x ~ s i)- 
Then 

(n n \ 

i=l i=l / 

Note that g alternates / and f ^ g (otherwise G = 0). Hence ^ILi s * < X^=i r «> which yields 
that c > since G is standard. It follows that b + d > from ad < be. Thus H = bf + dg is 
standard, and so is F = xG + H. The proof is complete. □ 

3 Applications and related topics 

In this section we apply the results obtained in the previous section to derive several known 
facts and to solve certain new problems in a unified manner. 

3.1 Orthogonal polynomials 

The classical orthogonal polynomials satisfy a three-term recurrence relation 

p n (x) = (a n X + b n )p n -i(x) - C n p n -2(x) 

with po(x) = 1 and p\(x) = a±x + b\, where a n ,c n > and b n £ R for all n > 1 (see |53t 
Theorem 3.2.1]). By Theorem 12.11 or Corollary 12.41 the orthogonal polynomials p n (x) form a 
Sturm sequence. This is a standard result in the theory of orthogonal polynomials. The following 
are some classical orthogonal polynomials (see Szego [53J for details). 

• (Tchebyshev) T n+ \(x) = 2xT n (x) — T n -\(x), T\(x) = x or T\(x) = 2x. 

• (Hermite) H n+ \(x) = 2xH n (x) — 2nH n -i(x), H\(x) = 2x. 

• (Laguerre) (n + l)L n+ i(x) = (2n + 1 — x)L n (x) — nL n -i(x), L\(x) = 1 — x. 

• (Legendre) (n + l)P n+ i(x) = (2n + l)xP n (x) - (n - l)P n _x(x), P\(x) = x. 

• (Gegenbauer) For A > —1/2, 

(n + l)C* +1 (x) = 2(n + X)xC^(x) - (n + 2A - l)C^(x), C?(x) = 2Xx. 
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• (Jacobi) Fora,/3> — 1, 

2n(n + a + /3)(2n + a + - 2)P^(x) 
= (2n + a + /3 - l)[(2n + a + /?)(2n + a + (3 - 2)x + a 2 - f3 2 ]P^\x) 
-2(n + a- l)(n + /3- l)(2n + a + /3)P^f } (x), n > 2, 

with P 1 (o '^ ) (z) = |(a + /? + 2)x + |(a - /?). 

For a = /3 = 0, -Pi ' "* (a;) reduces to the Legendre polynomials. The Gegenbauer polyno- 
mials and the Tchebyshev polynomials can also be viewed as special cases of the Jacobi 
polynomials. 

Note that the leading coefficients of the Laguerre polynomials L n (x) have the sign (— l) n . Set 
L n {%) = L n (—x). Then L n {x) are standard and satisfy 

(n + l)L n+ i(x) = (2n + 1 + x)L n {x) - nL n -x{x). 

Thus {L n (x)} forms a Sturm sequence, and so does {L n (x)}. 

3.2 Matching polynomials 

Let G be a graph with n vertices and p(G,k) the number of matchings of size k, i.e., the 
number of sets of k edges of G, no two edges having a common vertex. Set p(G, 0) = 1 for 
convenience. The matching polynomial M(G,x) = ^2 k (— l) k p(G, k)x n ~ 2k counts the matchings 
in the graph G. Clearly, for any v S V(G), 

P (G, k) = p(G - W, k) + £>(G - «}, k - 1), 

where the first term in the sum counts fe-matchings which do not use v and the second one 
counts /c-matching which do use v. This leads to a recurrence relation 

M(G, x) = xM(G - {v},x) - M (G - {v, u},x). 

From Theorem 12.31 and by induction on the number of vertices of graphs, it is easy to see 
that M(G, x) has only real zeros and M(G — {v},x) interlaces M(G,x) for any v E V(G), a 
well-known result (see Godsil and Gutman [25] for instance). 

The matching polynomials, formally introduced by Farrell [21] in 1979, have occurred not 
only in the combinatorial literature but also in various branches of physics and chemistry (see 
|25j for a brief survey). For example, to understand the behavior of a monomer-dimer system 
in statistical physics, Heilmann and Lieb [33] in 1972 introduced the partition function Q(G,x) 
of a monomer-dimer system, which is essentially the matching polynomial of the graph G with 
the edge weight W. If W = 1, then Q(G,x) is precisely the ordinary matching polynomial of 
the graph G. It is easy to yield the recurrence relation 

Q(G, x) = xQ(G - {v},x) - W ( u > V )Q(° ~ K 

Heilmann and Lieb \33\ Theorem 4.2] showed that if G is a simple graph with nonnegative edge 
weights, then Q(G, x) has only real zeros and Q(G — {v},x) interlaces Q(G, x) for any v G V(G). 
This result is now clear from the point of view of Theorem 12.31 
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The classical orthogonal polynomials are closely related to the matching polynomials. For 
example, the Tchebyshev polynomials of two kinds are the matching polynomials of paths and 
cycles respectively, the Hermite polynomials and the Laguerre polynomials are the matching 
polynomials of completes graphs and complete bipartite graphs respectively. 

Another class of polynomials closely related to the matching polynomials is the rook polyno- 
mials. The original rook polynomial defined in Riordan's book 00] counts the number of ways of 
arranging nonattacking rooks on a board (a board is a finite subset of N x N). Goldman, Joichi 
and White [26] conjectured that such rook polynomials have only real zeros. Wilf extended 
the concept of rook polynomial to the case of arbitrary matrix instead of a board (a board is 
identified with a (0, l)-matrix) and Nijenhuis showed that the extended rook polynomials of non- 
negative matrices have only real zeros by establishing a result similar to Theorem 12.31 (see [37] ) . 
Bender observed that the (extended) rook polynomial is the same as the matching polynomial 
of a (weighted) bipartite graph, and Nijenhuis's result and Goldman-Joichi- White's conjecture 
are therefore implied by the result of Heilmann and Lieb. 



3.3 Brenti's derangement polynomials 

Let S n denote the symmetric group on n-elements [n] = {1, 2, . . . , n}. Let tt be a permutation 
in S n . An element i € [n] is called an excedance of it if ir(i) > i. Denote by exc (n) the number 
of excedances of it. The permutation tt is called a derangement if 7r(i) ^ i for all i £ [n]. Let 
D n denote the set of all derangements of S n . Define the derangement polynomial 

dn(q)= E 9 CXCW - 

tt6-D„ 

For example, do(q) = l,di(q) = 0,d2(q) = q,dz(q) = q + q 2 - Since d n (l) = \D n \ we may consider 
d n (q) as a g-analogue of the derangement numbers. 

The exponential generating function of d n {q) can be written as 

-j-n -y 



n>0 



i-En>2(q + q 2 + --- + q n ~ 1 )t n /ri 



(see [TTJ HI] for instance). Gessel [21] gave a direct proof of (|3.ip with q = 1 based on a 
factorization of some so-called D-permutations. Kim and Zeng [35] gave a decomposition of 
derangements which interprets (|3.1|) directly. Using (|3.1|) Brenti [TT] showed that the polynomial 
d n {q) is symmetric and unimodal for n > 1. He further proposed the following. 

Conjecture 3.1 ( |11[ Conjecture]). The polynomial d n (q) has only real zeros for n > 1. 

The derangement polynomials are closely related to the classical Eulerian polynomials which 
are defined by 

Mq)= E g cxcW+1 

for n > 1 and Ao(q) = 1 (see Stanley |46t Proposition 1.3.12] for instance). Therefore, 



mq)/q = E E ? cxcW = E (£W«)- 

„r-r„i^cn, , i,— n N / 



sC[n] ttGZ?| s fc=0 

From the binomial inversion formula it follows that 



dn(q)=p(-lT~ k (fyA k (q)/q. 
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Recall that the Eulerian polynomials satisfy the recurrence relation 



A n (q) = nqAn-xiq) + g(l - g^.ife)' 



Hence {d n (q)} satisfies the recurrence relation 



d n {q) = (n - l)qd n -i(q) + q(l - q)d' n _ 1 {q) + (n - l)qd n - 2 {q). 



(3.2) 



Zhang [631 EU verified and generalized Conjecture 13.11 as follows: 

"Let f n {q) be a polynomial of degree n with nonnegative real coefficients and satisfy 

(a) f n (q) = a n qf n -i(q) + b n q(l + c^f^iq) + d n qf n _ 2 (q), where a n > 0, b n > 0, d n > and 
n > 2. 

(b) For n > 1, zero is a simple root of f n {q)- 

(c) fo(q) = e, = e\q and /2(g) has two real roots, where e > 0, e\ > 0. 

Then the polynomial has n distinct real roots, separated by the roots of /n-i (<?)•" 

However, Zhang's result does not hold in general and his proof is valid only for {d n (q)}. 
The condition c n < is necessary in Zhang's result, even in the case d n = 0. For example, let 
/o((Z) = hfl(q) = q, f2(q) = qfi(q) + q(q + 2)f{(q) = 2q(q + 1) and / 3 (g) = qf 2 (q) + q(2q + 
l)/^?) = 2q(5q 2 + 5q + 1) . Then zeros of f 2 {q) are —1 and and those of /3(g) are (— 5±\/5)/10 
and 0. Clearly, the latter cannot be separated by the former. 

Now from the point of view of Corollary |2.41 Coniecture l3.1l is obviously true and furthermore, 
the derangement polynomials d n (q) form a generalized Sturm sequence. 

Remark 3.2. The referee pointed out that Conjecture 13.11 has been proved by Canfield (unpub- 
lished). 

3.4 Narayana polynomials 

The Narayana numbers N(n,k) = -(£) ft."]) and the Narayana polynomials N n (q) = 
ELi iV ( n ' k )l k have many combinatorial interpretations and fascinating properties (see |50|. 
[5T| [52] for instance). For example, the Narayana numbers N(n, k) can be viewed as a refine- 
ment of the famous Catalan numbers C n since C n = N n (l) (see Remy [39] for a combinatorial 
proof of this and Stanley |47[ |4"9] for various combinatorial interpretations of the Catalan num- 
bers). Also, N n (2) are the Schroder numbers (see Foata and Zeilberger [23] for a combinatorial 
proof and Stanley [45] for an interesting history of the Schroder numbers). In fact, the original 
Narayana polynomials N n (q), introduced by Bonin, Shapiro and Simion [7], are defined as the 
g-analog of the Schroder numbers. It is shown in [7] that N n (q) has unimodal coefficients and 
has q = — 1 as one zero for n > 1 and q = — 2 as one zero for even n > 2. It is known that 



with N\(q) = q and N 2 (q) = q(l + q) (see [52] for instance). As pointed out by Stanley (see 
Bona [5]), Theorem 5.3.1 in Brenti [10] implies that N n (q) have only real zeros. Another proof 
for this result was recently found by Branden [9] by expressing N n (q) in terms of the Jacobi 
polynomials: 



N n {q) = N n (l + q) and 



(n + l)N n {q) 



(2n - 1)(1 + q)N n ^(q) - (n - 2)(1 - q) 2 N n ^ 2 



(q) 



(3.3) 
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More generally, Sulanke [50] defined polynomial sequences {p a ,n(q)}n>2, for <*. = 0,1,2, in 
terms of various "diagonal thickness" parameters on parallelogram polyominoes. These polyno- 
mials satisfy the recurrence relation 

(n + 1 - a)p atn +i(q) = (2n - 1 - a)(l + q)p a ,n(q) - (n - 2)(1 - g) 2 Pa,n-i(9) 

withp Qj 2( ( ?) = 9 andp ai 3(g) = g(l+g), and are generalizations of some well-known combinatorial 
sequences. In particular, po,n(q) = N n (q). 

Using Theorem 12.11 or Corollary 12.44 we can gi ve a niore direct and natural approach to the 
reality of zeros of N n (q) as well as p a ,n(q)- 

Proposition 3.3. Let {P n (x)} be a sequence of polynomials with nonnegative coefficients and 
degP n = degP n _i + 1. Suppose that 

P n (x) = (a n x + b n )P n -i{x) - (c n x + d n ) 2 P„_ 2 (x) 

where a n ,b n ,c n ,d n £ R. Then {P n (x)} forms a Sturm sequence. 

It is well known that the Catalan numbers (the Narayana numbers) count the number of 
1-stack sortable n-permutations (with k descents). The problem of stack sorting was introduced 
by Knuth [36] in 1960's and many variations have been considered since then (see Bona [6] 
for a survey). Let Wt(n,k) be the number of i-stack sortable n-permutations with k descents 
and W n f(q) = Sfc=o Wt( n > k)l k the associated generating function. In [5], Bona showed that 
{Wt(n, fc)}o<fc<n-i is a unimodal sequence for fixed n and t and further conjectured that W nt t(q) 
have only real zeros for n > 2. Note that W n i(q) = N n (q), the Narayana polynomial, and 
W n<n -i(q) = A n (q), the classical Eulerian polynomial. Hence the conjecture is true for t = 
l,n — 1. Branden [8] has recently verified the conjecture for t = 2,n — 2. But the conjecture 
remains open in the general case. 

3.5 Compositions of multisets and Dowling lattices 

Let n = (m, ri2, • • •) be the multiset consisting of rij copies of the ith type element. Denote 
by C(n, k) the number of compositions of n into exactly k parts. Then 

(nj + l)C(n + ej, k) = W(n, k - 1) + ( nj + k)0(n, k), (3.4) 

where n + e,- denotes the multiset obtained from n by adjoining one (additional) copy of the jth 
type element (see Riordan [40], p. 96]). Let f n (x) = X^fc>o ^( n 5 k) xk be the associated generating 
function. Using (|3.4p Simion [12] deduced the recurrence relation 

(rij + l)f n+ej (x) = (x + nj)f n (x) + x(x + l)fU x )- ( 3 - 5 ) 

By means of appropriate transformation to (|3.5p . Simion [4:2\ Theorem 1] can show that the 
polynomial f n (x) has all zeros in the interval [—1,0], and furthermore, f n (x) and / n + e .(x) have 
interlaced zeros. This result is now clear from the point of view of Theorem 12.11 

In particular, if n = (1,1,..., 1), then 0(n, k) = k\S{n,k) where S(n,k) is the Stirling 
number of the second kind. Thus the polynomial F n (x) = X^fc=i k\S(n, k)x k has only real 
zeros, which can also be followed from the formula F n (x) = ~^fT^ n (~r^) wn ere A n (x) is the 
Eulerian polynomial (see [5] for instance). The polynomial F n (x) was first studied by Tanny [54] . 
Benoumhani [2] gave a generalization of F n (x) replacing the Stiling numbers by the Whitney 
numbers of the Dowling lattices. The Dowling lattice Q n (G) is a geometric lattice of rank n over 
a finite group G of order m and have many remarkable properties (see [H H] [3] [20] for instance). 
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When m = 1, that is, G is the trivial group, Q n (G) is isomorphic to the lattice H n +i of partitions 
of an (n + l)-element set. So the Dowling lattices can be viewed as group-theoretic analogs of 
the partition lattices. Let W m (n, k) be the kth Whitney numbers of the second kind of Q n (G). 
Denote D m (n;x) = Xlfc=o Wm(n, k)x k and F n (m;x) = Ylk=o k\W m (n, k)x k . (Dowling gave a 
combinatorial interpretation for the coefficients k\W m (n,k), see [2].) Then 

D m (n; x) = (x + l)D m (n - 1; x) + mxD' m (n - 1; x) 

(see [3]) and 

-F n (m; x) = (x + l)F n _i(m; x) + x(x + m)F^__ l (m; x) 

(see [2]). Benoumhani showed that both D m (n;x) and F n (m;x) have only real zeros for n > 
1 (see [3j Theorem 2] and [2j Theorem 6] respectively). These results can also be followed 
from Corollary 11.21 As a consequence, W m (n, k) is unimodal and log-concave in k. This gives 
supports to a long-standing conjecture that the Whitney numbers of the second kind of any 
finite geometric lattice are unimodal or even log-concave (see [44, Conjecture 3]). When m = 1, 
we have Q n (G) = H n +i and W m (n, k) = S(n, k). Again we obtain that the polynomials S n {x) = 
J2k=o ^( n ' k)x k and F n (x) = J2k=o k\S(n, k)x k have only real zeros for n > 1. 

3.6 Eulerian polynomials of Coxeter groups 

Given a finite Coxeter group W, define the Eulerian polynomials of W by 

P(W,x) = x dw{n \ 
new 

where diy( 7r ) is the number of VF-descents of tt. We refer the reader to Bjorner and Brenti [1] 
for relevant definitions. Brenti proposed the following conjecture. 

Conjecture 3.4 ( |13l Conjecture 5.2]). For every finite Coxeter group W , the polynomial 
P(W; x) has only real zeros. 

For Coxeter groups of type A n , it is known that P(A n ,x) = A n (x)/x, the shifted Eulerian 
polynomial. The classical Eulerian polynomial A n (x) has only real zeros, so does P(A n ,x). 
Since {A n (x)} satisfies 

A n (x) = nxA n _i(x) + x(l - x)A' n _ 1 (x), A (x) = 1, 

it immediately yields that {P(A n ,x)} satisfies 

P(A n ,x) = (nx + l)P{A n ^ x ,x) + x(l - x)P'(A n _ l ,x), P(A , x) = 1. 

In [22], Foata and Schiitzenberger introduced a g-analog of the classical Eulerian polynomials 
defined by 

A n {x;q) = ^ 

vr€<S„ 

where exc (tt) and c(ir) denote the numbers of excedances and cycles in ir respectively. It is 
clear that A n (x; 1) = A n (x) is precisely the classical Eulerian polynomial. Brenti showed that 
g-Eulerian polynomials satisfy the recurrence relation 

d 

A n (x;q) = (nx + q- l),4 n _i(x; 9) + x(l - x)^-A n _ 1 (x;q), 
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with Aq(x; q) = x ([16] Proposition 7.2]). He showed also that A n (x; q) has only real nonnegative 
simple zeros when q is a positive rational number Theorem 7.5]). 

For Coxeter groups of type B n , Brenti [T3] defined a g-analogues of P(B n ,x), which reduces 
to A n {x) when q = and to P(B n , q) when q = 1, by 



where N(tt) = \{i £ [n] : ir(i) < 0}|. He showed that {B n (x; q)} satisfies the recurrence relation 



with Bo(x;q) = 1 ( |13l Theorem 3.4 (i)]) and that all B n {x;q) have only real zeros for q > 
( \13\ Corollary 3.7]). In particular P{B n ,x) has only real zeros. 

By the classification of finite irreducible Coxeter groups, it suffices to decide whether the 
conjecture holds for Coxeter groups of type D n to settle Conjecture 13.41 (see Brenti [13] for 
further information). 

Using Theorem 12.11 or Corollary 12.41 we can give a unified interpretation of the reality of 
zeros of A n {x), P(A n , x), A n (x; q) and B n (x;q). 

Proposition 3.5. Let {P n (x)} be a sequence of polynomials with nonnegative coefficients and 
degP n = degP n _i + 1. Suppose that 



where a n , b n E M and c n < 0, d n > 0. Then {P n (x)} forms a generalized Sturm sequence. 

3.7 Genus polynomials of graphs 

Given a finite graph G and a nonnegative integer k, let ^y(G,k) denote the number of 
distinct embeddings of the graph G into an oriented surface of genus k. We refer the reader to 
|29| for the basic terminology of graph embeddings. The genus polynomial is defined in [27] as 
GP(G, x) = X^fc>o t(^) k)x k . Gross, Robbins and Tucker [28] showed that the genus distribution 
of the bouquet is log-concave and conjectured that the genus distribution of every graph is log- 
concave (i.e., j(G,k) is log-concave in k). Stahl [43], Conjecture 6.4] further conjectured that 
the genus polynomial of every graph has only real zeros. He also verified the conjecture for 
several infinite families of graphs by establishing the recurrence relations of the associated genus 
polynomials. These results follow from our results in the previous section. In particular, Stahl 
considered the ff-linear family of graphs obtained by consistently amalgamating additional 
copies of a graph H. For such a family {G n }, there is a square matrix M and a vector v with 
entries in Z[x] such that the genus polynomial of G n is the first entry of M n v ( |431 Proposition 
5.2]). The following are genus generating matrices and initial vectors of certain linear families 
given by Stahl [4*3] . 

Example 3.6. 



B n {x;q) = ^ W ^ flW 



7TGB; 



B n (x; q) = {l+ [(1 + q)n - l]x}S n _i(x; q) + (1 + q)x(l - x)—B' n _ 1 (x; q) 



P n (x) = (a n x + b n )P n -i(x) + x(c n x + d n )P^_ 1 (x) 
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/ 3x 3 \ ( 2 

(iii) For the double ladders, M3 = 6 I and = 2 



2x l + 3xy J V 1 + a; '' 

(iv) For the diamonds, M4 = 4 ( ^ ^~ ^ J and = 2 ( ^ X 

^ t? ^ + • ! ! ^ Af f 192x 96 + 288x\ /l8 + 18x\ 

(v) For the triple ladders, M 5 = + ^ ^ + ^ j and v 5 = ^ fl + 3Qx j 

8 + 68x 4 + 16x\ / 2 + 14s 

32x + 48x 2 16x J V6 " 1 8s + 8x 2 



(vi) For the ^-linear graphs, M 6 = ^2 i«„ ) and ^6 

(vii) For the W4-linear graphs, 

_ (2 + 65x + 54x 2 1 + 22x \ /2 + 58s + 36s 2 

Ai7 " 4 V 16s + 104x 2 8x + 16x 2 ) and U7 " V 16 + 80s 

(vih) For the triangular prisms, 

/ 162x 54 \ / 8 \ 

M 8 = 24x 2 72x 12 + 108x and v 8 = 4 + Ax . 

\lls 2 15x + 117s 2 l + 72s / \l + 7s/ 

Stahl verified his conjecture for the linear families associated with Example 13.61 (i), (ii) and 
(iv), i.e., the first entry of M k v has only real zeros for k > 1. He left the remaining as a 
conjecture and asked some more general questions as follows. 

Conjecture 3.7 ([43j Conjecture 6.9]). The zeros of the genus polynomials of the graphs listed 
in Example \3.6\ are real and negative. 

Question 3.8 ( |43^ Question 6.10 and 6.11]). Let M(x) be a square matrix whose entries are 
real polynomials. 

(i) Under what conditions, if (f(x),g(x)) is a pair of polynomials whose zeros interlace, do 
the zeros of the two components of the vector (f(x),g(x))M(x) interlace? 

(ii) Under what conditions, are the zeros of each of the entries of M k (x) all real for k = 
1,2,...? 

(2 + 3s 1 \ 

Stahl showed that the matrix M4 = 4 ( J in Example 13.61 has both properties in 

/% x 3 \ 

Question 13.81 but the matrix M = I I has neither. In what follows, we apply 

\ ZiX oOC ~\~ 1UUU J 

Theorem 12.61 to give an answer to Question 13.81 and to verify Conjecture 13.71 for the graphs in 
Example 13.61 (i)-(vi) in a unified approach. 

(a(x 1 c(x 1 \ 
. , I , { J , where a, b, c, d are polynomials with nonnegative 

coefficients. We say that the polynomial matrix M is nice if 

(a) dega,deg<i < l,deg6 < 2 and c is a positive constant. 

(b) det(M) > for x < 0. 
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Proposition 3.10. Let M be a nice matrix. Suppose that f,g € PF and g ^ j '. Then the 
following hold. 

(i) If (F, G) = (f,g)M, then F, G £ PF and G<F. 

(ii) //(Gi,F;l) t = M(g,f f, then F 1 ,G 1 G PF and Gi ■< F 1 . 
(Hi) Each entry of M k has only real zeros for k = 1, 2, . . .. 

Proof. We have F = af + bg and G = cf + dg. We may assume that F ^ and G ^ 0. To 
prove (i) by means of Theorem 12.61 it suffices to prove that degF = degG or degG + 1. This is 
obvious if deg/ = degg + 1. So let deg/ = degg = n. We distinguish two cases. Suppose first 
that degG = n + 1. Then degd = 1. Since ad — be > whenever x < 0, we have deg a > 1 or 
deg b > 1, Hence degF > n + 1. On the other hand, it is clear that degF < n + 2 since deg a < 1 
and deg 6 < 2. Suppose now that degG = n. Then d is a constant. Again by the assumption 
ad — be > whenever x < 0, we have deg b < 1. It follows that n < degF < n + 1. Thus (i) 
follows from Theorem 12.61 

Similarly, (ii) follows from Theorem 12.61 since F\ = df + bg and G\ = cf + ag. 

Finally, we apply (i) to prove (iii). For k = 1, 2, . . ., let 

Then (Z^,^) = (a,c),(fi 2 \g { ^) = (b,d) and (f^g^) = (/f )M for i = 1,2. By the 
assumption ad — be > whenever x < it follows that d(r) = implies b(r) < 0. This means 
that b £ PF and d ^ b. Also, c ^ a since degc = and deg a < 1. Thus (iii) follows from (i) by 
induction on k. □ 

Proposition 3.11. The zeros of the genus polynomials of the graphs listed in Example \3.6\ (i)- 
(vi) are real and negative. 

Proof. We need to verify that for i = 1,2,... ,6, the first entry of the column vector M^Vi has 
only real zeros for each k. By Proposition 13.101 (ii) and by induction on k, it suffices to prove 
that each Mi is the product of certain nice matrices and ^ where v% = {v^\vf^) T . 

For i = 1, 2, 3, 4, it is easy to verify that each matrix Mj is nice. For i = 5,6, the matrix 
can be decomposed into the product of certain nice matrices: 



M 5 = 24 

Ma = 



4 + I2x 8 W0 1 
1 + 12a; 3 + 8x) \x 

1\ /8 + 12x 4 
x oj V2 + 17x l + 4x 

It is also easy to verify that ^ for i = 1, 2, . . . , 6. Thus the proof is complete. □ 

Remark 3.12. The genus polynomials of the W^-linear graphs in Example 13.61 (vii) may have 
nonreal zeros. For example, let uj = Mjvj. Then the first entry of u-j is 



10 + 339x + 2855x 2 + 2736x 3 + 972x 4 ) 



Using Mathematica, we can obtain the approximations of four zeros of : 

x^l = -1.34194 ± i0.88376, = -0.0828403, xf> = -0.0481022. 
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This gives a counterexample to Conjecture 13.71 Stahl's conjecture about the reality of zeros of 
the genus polynomials is therefore false in general. But it is possible that the genus distribution 
of each graph is log-concave. 

We end this paper by proposing the following. 

Problem 3.13. Characterize all real polynomial matrices that can be decomposed into the prod- 
uct of finite nice matrices and find an algorithm of decomposition for such matrices. 
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